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Abstract 

The AdS/hydrodynamics correspondence provides a 1-1 map between large wavelength 
features of AdS black branes and conformal fluid flows. In this note we consider boundaries 
between nonrelativistic flows, applying the usual boundary conditions for viscous fluids. 
We find that a naive application of the correspondence to these boundaries yields a surface 
layer in the gravity theory whose stress tensor is not equal to that given by the Israel 
matching conditions. In particular, while neither stress tensor satisfies the null energy 
condition and both have nonvanishing momentum, only Israel's tensor has stress. The 
disagreement arises entirely from corrections to the metric due to multiple derivatives of 
the flow velocity, which violate Israel's finiteness assumption in the thin wall limit. 
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1 Introduction 



It has long been known that the dynamics of a p- dimensional gravitational theory is 
captured by quantities on (p — l)-dimensional hypersurfaces [JJ. It was argued by Damour 
[2], based on an analogy by Hartle and Hawking [3J, that in the case of certain black hole 
solutions these surface quantities describe the flow of a viscous (p — l)-dimensional fluid. 
Perhaps the most concrete realization of this idea is the one to one map between large 
wavelength features of asymptotically AdS p black brane solutions and (p— l)-dimensional 
conformal fluid flows presented recently in Refs. [I] and [5]. 

This AdS /hydrodynamics correspondence provides an explicit black brane solution for 
every history of a particular conformal fluid, so long as the fluid variables are constant 
over distances large compared with the inverse temperature. For example progress towards 
gravity duals of shock waves and vortices has appeared in Refs. [6] and [7j. In particular, 
there must be gravity duals to turbulent flows. Turbulence is generic in fluid flows under a 
wide range of conditions. The dual of these fluid conditions then provides some condition 
on a gravity solution under which it to generically decays into a turbulent configuration. 
An example of such a situation was presented in Ref . [8] . It would of course be interesting 
to characterize the gravity duals of turbulent flows, and of the conditions under which 
turbulence may be expected. In hydrodynamics, even the most basic scaling laws are 
altered by turbulence. If gravitational solutions near, for example, spacelike singularities 
(where indeed chaotic evolution is expected [9]) or certain event horizons do generically 
decay to turbulent solutions, it would be difficult to overstate the potential consequences 
for, for example, the horizon problem. 

Perhaps the best understood turbulence is steady state turbulence, in which energy 
is injected into a system at the same rate at which it dissipates. Richardson's cascade 
model [10] of steady (3+l)-dimensional turbulence is as follows. Energy is injected into 
a system at large characteristic distance scales, for example, a lake warms the air. This 
creates large vortices, which decay into smaller vortices. Thus the energy flows to smaller 
distance scales. At sufficiently small distance scales, higher order derivative terms in the 
equations of motion become relevant, such as viscosity terms. These lead to dissipation 
of the energy in sufficiently small vortices. Thus energy cascades from the long length 
scale in which it is introduced, down to the dissipation scale. 

To realize steady state turbulence, one needs to inject energy into a system. There are 



1 



two principal ways to do this. First, one may deform the fluid via external perturbations. 
Second, one may apply boundary conditions, for example one may consider fluid flow 
in a pipe or wind tunnel. The first approach was applied to the AdS/hydrodynamics 
correspondence in Ref. [8], where it was argued that a laminar fluid flow and the dual 
gravity solution decay to turbulent configurations. This approach has the disadvantage 
that solutions are quite complicated, due to the necessarily inhomogeneous forcing and 
to the geometric implementation of the forcing on the gravity side. 

In this note we will take a preliminary step towards a realization of the second ap- 
proach to creating steady state turbulence, we will investigate boundary conditions in 
the AdS/hydrodynamics correspondence. For simplicity, we will consider nonrelativistic, 
incompressible flows. Consider the surface which separates a solid object from such a 
fluid. The normal velocity of the fluid into the solid must vanish. If furthermore the fluid 
is viscous, as fluids in the AdS/hydrodynamics correspondence are [2], then the tangential 
relative velocity of the fluid must also vanish. 

What does this correspond to on the gravity side? The answer to this question is not 
necessarily unique, one may define a dual and then attempt to understand its dynamics. 
One interesting case, which is already sufficient to generate turbulence, is a solid which 
is a thin, infinite sheet with a stationary fluid on the left side and a moving fluid on the 
right. In this case a natural choice would be to consider the gravity duals of both fluids 
and then to attempt to glue them together. Equivalently one may choose to think of the 
entirety of the left side as a solid wall, filling the left half of spacetime, and a liquid filling 
the right half. The wall is stationary and so one chooses the dual to be a stationary black 
brane in half of AdS. Whatever one chooses to think, the logic is that one imposes that 
the left half of the gravity dual be a static black brane in AdS, and that the right side be 
the gravity dual given by the prescription of Ref. [I] . 

So how does one glue these two vacuum gravity solutions together? Clearly there are 
many inequivalent choices. One possibility is to simply attach them and then use the 
Israel matching conditions [llj to determine the stress tensor on the surface layer that 
separates the two sides. This is equivalent to letting the gravitational solution continu- 
ously interpolate between the two solutions over a finite distance d and then taking the 
limit as this distance tends to zero. While there are many ways of performing this inter- 
polation, so long as the extrinsic curvature is kept finite, they all lead to the same stress 
tensor as the interpolation distance d — > 0. 

Another possibility is to let the fluid configuration continuously interpolate between 
the two solutions, and then take the dual using the prescription of Ref. [I] . As the fluid is 
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not a solution of the Navier-Stokes equation in this region, the dual will not be a solution 
of the vacuum Einstein equations in this region. Instead it will solve Einstein's equations 
with a nonvanishing stress tensor supported on a surface layer. The ultralocality of the 
duality map implies that the vacuum Einstein equations will however be solved away 
from the surface layer. In this case, one cannot take the interpolation distance d to zero, 
because the dual is not defined when derivatives are large with respect to the inverse of 
the temperature T. Thus the minimum size of d will be of order 1/T. Again there are 
many inequivalent ways of performing the interpolation. But we will see that, at least 
for the quantities at we are able to calculate, when d is large with respect to 1/T, the 
difference between these prescriptions is suppressed by powers of dT and so, like Israel's 
method, there is a single answer. 

The perhaps surprising result is that the two methods yield bulk stress tensors which 
differ by a finite amount. They did not need to agree, indeed one is derived at small d and 
the other for large d. The reason that they disagree is as follows. The construction of the 
metric from the fluid flow proceeds order by order in the derivatives of the fluid's velocity. 
The boundary conditions imply that the velocity of the fluid is the same on both sides of 
the wall, however the first derivatives differ. Therefore, whatever regularization scheme 
one uses on the fluid side, the second derivative of the velocity diverges at small d. This 
means that the metric corrections derived using the map of |4J will diverge at small d, 
invalidating the finiteness assumption in Israel's derivation. In fact, we will see that the 
disagreement between the two calculations of the stress tensor differ only in these higher 
derivative terms. Of course the fluid map is not defined at small d, as it yields a divergent 
series, and so no divergences appear within the range of validity of either approach. 

We will begin in Sec. El by describing the flow of interest. The velocity will be kept 
sufficiently arbitrary to allow a general interpolation between the flows on the two sides 
of the wall, and in Sec. El the naive gravity dual will be calculated using the prescription 
of |Q. We will see that those higher order derivative corrections which we calculate are 
indeed suppressed by factors of dT. Then in Sec. H] we will calculate the bulk stress 
tensor of the interpolation between the two gravity solutions. First it will be calculated 
for the interpolation dual to a continuously interpolating fluid flow. It will be seen that 
contributions from the second derivative of the velocity are ^-independent, while higher 
order contributions are suppressed by powers of dT. Thus the result is independent of the 
interpolation scheme when d is sufficiently large. The stress tensor will then be calculated 
directly from the Israel matching conditions on the two solutions of the vacuum Einstein 
equations. It will be seen that the two stress tensors agree up to terms corresponding to 
a divergence in the extrinsic curvature at small d, and that only the second stress tensor 
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contains a nonvanishing stress. 



2 The Flow 



2.1 The ansatz 

We will consider a hydro dynamic flow in 4-dimensional Minkowski space, using a (— , +, +, +) 
metric. To highlight the essential features of the boundary condition, we will consider 
the simplest possible flow. The liquid will only move in the y direction, with a velocity 
v = v(x) that only depends on the coordinate x. The velocity will be taken to be small, 
and we will drop all terms which are quadratic in v. In fact, as described in Refs. [121 E] 
we will work in the nonrelativistic, incompressible limit. More precisely, we will show that 
our flow satisfies both the full relativistic equations of motion at order 0(v) and also the 
incompressible Navier-Stokes equation. 

We will set c = 1 . The conformal fluid which is dual to Einstein gravity with a negative 
cosmological constant is very particular. Being conformal, all of its transport coefficients 
may be expressed in terms of a single dimensionful quantity, such as the temperature T, 
and certain constants which may be calculated from the gravity dual. In the case at hand 
for example the shear viscosity rj, pressure p and density p have been found in Ref. [4] 

^i4 t4 - p ^ Ti (2 - i) 

where Gat is the dual Newton's constant. 

The relativistic velocity 4-vector u is, to linear order in v, simply 

<V = ( y=L= , 0, 7 =L= , 0) ~ (1, 0, v, 0). (2.2) 



vr^ 1 vt 



V 



We will be interested in the fluid velocity in three regions, as illustrated in Fig. [TJ First, 
on the left, where v = 0. Second, we will be interested in the velocity on the right, where v 
will be linear in x. We will show momentarily that this is a solution to the hydrodynamic 
equations of motion and so will be dual to a vacuum solution of Einstein's equations. 
Finally, we will be interested in an interpolating region where v will be arbitrary and 
we will not impose the equations of motion, therefore the dual metric will not solve the 
vacuum Einstein equations but, like any metric, will solve Einstein's equations with some 
stress tensor. 

Clearly the left, v — 0, satisfies the fluid equations of motion. We will now verify that 
the region on the right satisfies the relativistic equations of motion, which are simply the 
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Figure 1: The fluid velocity v is in the y direction, and it depends on the x coordinate. 
On the left the fluid is stationary, on the right the fluid velocity is linear. These two 
regions solve the fluid equations of motion at linear order in v. There is an interpolating 
region of width d, which must be larger than the inverse temperature, in which v does 
not satisfy the equations of motion, v and its first derivative v' are continuous at x = 
and x = d. 

conservation of the stress tensor 



In accordance with the usual fluid approximation [13J, we will work at large enough dis- 
tance scales that only the velocity v and its first derivative v' need be considered in the 
stress tensor. This approximation in general is problematic, leading for example to su- 
perluminal propagation [TJ]. However, as we will be interested in velocities well below 
the speed of light, no problems will arise. Later, when we will consider the interpolating 
region, where the second derivative may be large, we will make no such approximation. 
We will consider the bulk stress tensor to higher order, calculating all terms up to two 
derivatives and several terms up to three or four derivatives to check that they are sub- 
dominant. However we do not impose that the interpolating region satisfies the equations 
of motion, indeed that would lead to a vanishing bulk stress tensor. 

2.2 Relativist ic and nonrelativistic equations of motion 

Dropping all higher derivatives of the velocity and using the fact that the fluid is conformal 
to eliminate the bulk viscosity and replace p with 3p, the hydrodynamic stress tensor is 




(2.3) 



T» v = v {rf v + Au»u v ) - 2r}<7> 



(2.4) 
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where the shear strain rate a^ v is defined as 



-d x u x P^. 



(2.5) 



Here parenthesis denote symmetrization with a factor of one half and P^ u is a projector 
onto the spacelike directions in the reference frame of the fluid 
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(2.6) 



Substituting the velocity ansatz ( 12. 2 p into the definition (I2.5P one easily finds the shear 
strain at linear order in v 
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(2.7) 



The constants of proportionality ( 12.1 p in this particular fluid can then be inserted into 
the general formula ( 12 .4p for T^ w to express the stress tensor in terms of the temperature 
T and the velocity v 
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The velocity only depends on the coordinate x. Let us choose boundary conditions so 
that the temperature T also only depends on x. Then the equations of motion (I2.3P are 
simply 

= d x T xv . (2.9) 

However all of the components T xu are constants except for T xx and T xy . Therefore the 
only nontrivial equation of motion at linear order in v is 



= d x T xx + d x T xy 



7T 3 T 3 



T' 



71 



-d x (T 3 v'). 



(2.10) 



AG N 16G N 

When v is linear in x, its second derivative vanishes. Therefore this equation of motion, the 
conservation of momentum in the y direction, implies that the temperature is constant. 
More precisely it implies that (dT)/T is negligible in this approximation. In light of 



6 



the relations ( 12. ip this is consistent with the incompressibility assumption that we have 
imposed on our fluid. 

Therefore we recover the fact that at sufficiently small velocities, a constant tempera- 
ture and a y- velocity which is linear in x solve the equations of motion (12.31) . Intuitively 
this is clear. Further to the right, the fluid is moving faster. Therefore the viscous force on 
a unit of fluid exerted by the faster fluid on its right (in the +x direction) will accelerate 
it in the +y direction, whereas the slower fluid on its left will exert a viscous force that 
decelerates it. A steady flow occurs when these two forces cancel, which implies that the 
second derivative of the flow vanishes. Had there been a temperature gradient, then the 
viscosity to density ratio would also have been stronger on one side by (12. ip . and so this 
balance could only be maintained by introducing a second derivative of the velocity. 

Clearly a linear velocity also satisfies the nonrelativistic Navier-Stokes equation for an 
incompressible, Newtonian fluid 

p(d t v k + (v ■ d)v k ) = -d kP + uV 2 v k . (2.11) 

In fact, each term vanishes independently. Note that the vanishing of the dp term is not 
merely a consequence of incompressibility. In incompressible flows it may be of the same 
order as the viscous term. It vanishes in this case because this provides a solution to ( 12.111) 
and it is consistent with the various nonrelativistic, small gradient and incompressible 
limits taken above. 

In conclusion, we have considered fluid flows with an x-dependent velocity v in the y 
direction. We have verified that, to linear order in v, these satisfy both the relativistic 
and nonrelativistic equations of motion when v is linear in x and the temperature T is 
constant. Our flows of interest will have v = on the left, v linear on the right and an 
interpolating region inbetween. Thus the equations of motion will be satisfied on the left 
and the right but not in the interpolating region, leading to a dual gravity solution which 
solves Einstein's vacuum equations on the left and right, but inbetween requires material 
described by a nontrivial stress tensor. We have also found formula (12.61) and (12. 7p for 
the projector and the shear tensor a^ v for general functions v, and so these results 
may be applied to the interpolating region. 

Clearly if the fluid velocity is linear over a large enough distance, it will eventually 
approach the speed of light and the nonrelativistic approximation will break down. There- 
fore our analysis is only relevant near the boundary. The solution may be made global 
by introducing a second boundary, such that the velocity is constant on the other side of 
the second boundary. We will see below that the stress tensor on the second boundary, to 
linear order in v, will be minus the stress tensor of the first boundary. Of course at higher 
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orders one may expect an attraction between the two boundaries, and so this solution will 
not be stationary. The underlying assumption in this note is that the walls on the gravity 
side are built of a solid which remains fixed. The fact that the null energy condition is 
violated may be a sign that such a material would be inconsistent, as some null energy 
violating configurations lead to superluminal propagation or instabilities [15], although 
some do not [TB], in which case this configuration should only be considered over a suffi- 
ciently short timescale. It may be that this timescale is never sufficient for turbulence to 
develop. 

3 The Gravity Dual 

The AdS/hydrodynamics correspondence yields a black brane metric dual to arbitrary 
flows in very particular conformal fluids, which for example obey the relations (12. ip . If the 
flow satisfies the hydrodynamic equations of motion (12. 3p . the dual satisfies the vacuum 
Einstein equations, in this case with a negative cosmological constant. 

3.1 A note on ultralocality 

The correspondence, at least in the incarnation in Refs. [I] and [5], is ultralocal. This 
means the following. Consider the set of ingoing null geodesies which run from the bound- 
ary to the black brane horizon. Clearly each point in the bulk is on precisely one such 
geodesic. Also each point on the boundary is on one such geodesic. Therefore these 
geodesies can be used to associate a fixed single boundary point to each bulk point. This 
association is not one to one, there is an entire geodesic worth of bulk points associated 
to each boundary point. 

To implement the map, one identifies the boundary with the Minkowski space on 
which the fluid lives. The metric and its derivatives at a point in the bulk are determined 
entirely by the fluid velocity and temperature and their derivatives at the associated 
boundary point. One does not need to know the behavior of the fluid elsewhere. This is 
the ultralocality of the correspondence. In particular, the fact that the fluid satisfies the 
hydrodynamic equations of motion on the left and right (at x < and x > d) implies that 
the dual metric will satisfy the vacuum Einstein equations on the left and right, so long 
as the characteristic distance over which these quantities vary is greater than the inverse 
temperature. 

Our fluid does not satisfy the hydrodynamic equations of motion at < x < d. 
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This means that the dual gravitational configuration will not satisfy the vacuum Einstein 
equations, instead it will only satisfy Einstein's equations with a nonzero stress tensor, 
which we will calculate in Sec. HJ 

This ultralocality is somewhat different from the ultralocality that one encounters in 
classical field theories, or in the BKL limit of gravity theories, in that the bulk geometry 
is ultralocal in terms of null and not temporal evolution. This ingoing null identification 
identifies the temporal evolution of the fluid with outward radial evolution for the gravity 
theory. That is to say, the metric at larger radii but the same time is determined by the 
fluid in the future but at the same location. In particular, a timeslice of the bulk geometry 
is determined by the evolution of the boundary fluid during a fixed interval of time. This 
interval is of order the inverse temperature, and so no appreciable evolution may occur 
during this interval if the temperature is large enough for the correspondence to hold. In 
this sense any fixed timeslice of the gravity dual contains only as much information as a 
fixed timeslice of the fluid, despite being one dimension greater. 

As each event in the fluid is identified with an inward null geodesic in the bulk, the 
metric corresponding to this event appears to be falling towards the black hole at the 
speed of light. This is not at all to say that there is a Killing vector in the inward null 
direction, the metric changes in that direction, but in a fashion which is fixed by the map. 
Thus a disturbance on the boundary creates gravity waves which fly inward at the speed 
of light to the horizon. Similarly, pasting together an infinite sequence of bulk timeslices 
which are separated by time intervals 1/T, one obtains a pattern which falls from the 
boundary in to the horizon at the speed of light. Although each individual timeslice is 
too small to see any evolution, the entire pattern is dual to the entire history of the flow. 
Like a movie reel, the pattern in turn allows one to reconstruct the gravity dual, as it 
contains the timeslices. 

One may use this identification to speculate on the gravitational dual of decaying 
turbulence. For example, the inverse cascade of decaying (2+l)-dimensional turbulence 
consists of a chaotic period during which the fluid is subjected to random external forces 
followed by a relaxation period, characterized by the merging of well-separated vortices 
[T7J [18]. This would then be dual to a kind of forest of gravity waves falling from the 
boundary to the horizon, beginning when the boundary is subjected to a random pertur- 
bation. First the canopy, representing the chaotic period, falls out of the boundary into 
the horizon. When the external perturbation is turned off, it is followed by the branches 
representing the vortices, and then the branches merge into trunks as the vortices merge. 
The usual cascade [TH], [2U] in (3+l)-dimensional turbulence may be similarly described, 
but when the boundary is randomly perturbed, the trunk falls out first. This leads to 
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the rather bizarre observation that black brane geometries in A0IS4 and AdSs respond 
very differently to random perturbations of their boundaries. Needless to say, it would 
be interesting to make this picture precise, or to see whether it is inconsistent with the 
various approximations involved in the duality. 

3.2 The metric 

We will now calculate the metric dual to the flow (12. 2p using the map in Ref. jl] with the 
simplified notation of Refs. [5j |2T] . This map takes regions in which the flow satisfies the 
fluid equations to regions in which the metric satisfies the source-free Einstein equations. 
Acting on the region in which the fluid does not satisfy the hydrodynamic equation, the 
map is not known to have any special properties other than continuity, which will produce 
an interpolation between the vacuum Einstein metrics on the two sides. Therefore the 
choice of this map corresponds to a rather arbitrary choice of interpolation. However 
we will see that this interpolation has two nice properties. First, it is reasonably inde- 
pendent of the interpolating velocity function chosen. In particular, the third and higher 
derivatives of the velocity will yield contributions to the integrated stress tensor which are 
suppressed by powers of dT, while the leading contribution is independent of d. Second, 
the resulting stress tensor is simpler than the Israel stress tensor, it will have zero stress, 
whereas Israel's stress tensor has shear stress. 

Ultralocality in the ingoing null direction implies that the simplest coordinates in 
which to express the metric are Gaussian null coordinates, in which r parametrizes the 
ingoing null lines. In these coordinates, the bulk 5-dimensional metric corresponding to 
an x- dependent 4- dimensional fluid flow is [5] 

d s 2 = G MN dX M dX N = -2u„(x) dx» (dr + V u (r, x) dx u ) + &„ u (r, x)dx^dx u (3.1) 

where, up to second derivatives in v, V u and are defined as 



A 



1 (6r)P;P A a ! 



A 



v 




(3.2) 



and 



r 2 P llu -uluj Xu + 2{br) 2 Q l {br) 



<V + Qi(br)a a x , 



.1/ 



02(H 




d- 1 



(3.3) 



10 



The functions that appear in this definition are defined in Ref. [5], we will quote the 
definitions as they are needed. 

Eq. ( 13. ip is the bulk metric dual to a fluid flow in an arbitrary curved space. The 
fluid is conformally invariant, and the conformal invariance has been used to write the 
metric in a compact form using objects which transform covariantly under the conformal 
symmetry. We are interested in a flat boundary, and so many of these objects will vanish. 
In fact, since u only depends on x, but only has nonvanishing components in the t and y 
directions, even the gauge field for a Weyl transformation will vanish 

1 



d-1 







(3.4) 



as u A V xu^ = and V x u x = 0. This implies that the Weyl-covariant derivative reduces 
to the ordinary derivative 
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The Weyl-covariant Schouten tensor S is proportional to the Weyl-covariant curvature 
of the boundary. As the Weyl-covariant derivative is just the ordinary derivative, this is 
just the ordinary curvature. As the boundary is Minkowski space, the curvature vanishes, 
and so the Weyl-covariant Schouten tensor also vanishes 
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Similarly the Weyl-covariant Weyl curvature C is the sum of the ordinary Weyl curvature 
and the curvature of the Weyl tensor, which both vanish and so 
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The vorticity u does not vanish, however like the shear strain a it is of first order in 
v. Therefore u 2 , ua and a 2 terms are all of order 0(v 2 ) and so do not contribute at order 
0(v). Thus only the third and fourth terms of (13. 2\i contribute to V M . The third term is 
easily evaluated 
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The functions T are easily expressed in terms of the shorthand notation <>, which 
symmetrizes and contracts with the projectors P^ v . The first two are identically zero, 
while the second two are of order 0{v 2 ) 



= 2u a T> a a IMU = 
X 2At „ = C^pifvP = 

T 4 ^ = 4*7°^ ~ 0. 
Therefore only the first and third terms of (13.31) contribute to 0^ 
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ds 



Finally, inserting Eqs. (13 .9p and (I3.14p into (13. ip . one finds the final form of the metric 

-r 2 (l - -L) dt 2 - Idtdr + (2^ + ^(br)v") dtdy 

+ r 2 (dx 2 + dy 2 + dz 2 ) + 2r 2 6 Qi(br)v' ' dxdy - 2vdydr (3.15) 



where b = 1/ttT and the functions t>i and $ji are defined in Eqs. (14. 5 p and (14.61) . Using 
the basis (t, x, y, z, r), we may write the metric in matrix form as 
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3.3 Christoffel symbols 



In Sec. H] we will see that the leading contribution to the stress tensor comes from the 
second derivative of the velocity. Contributions at that order come from the second 
derivative of the velocity in the curvature, which in turn contains contributions from first 
and second derivatives of the velocity in the Christoffel symbols, as well as from Christoffel 
symbols which are velocity independent as these are multiplied by velocity-dependent 
terms when calculating the curvature. We will now calculate all of the Christoffel symbols 
up to first order in v, v' and v", although the v terms will not contribute to the stress 
tensor. 

We will begin with the terms at order 0(v°), these are just the Christoffel symbols of 
the static black brane 
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The new terms, at order 0(v), are 
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Ky = r r yx = ~\r 2 V - (r 2 - -L) (r& fll (&r)t/ + \rW &(br)j) 

= K* = l W + 2rb Ql (br)v> + rV Q [(br)v') . (3.17) 

3.4 The Riemann tensor and the Ricci tensor and scalar 

Using the Christoffel symbols one can now easily compute the Riemann tensor. The order 
O(v ) terms again are just those of the static AdS black brane 

3 

R trtr = l ~^i ( 3 - 18 ) 

Rtxtx = Rtyty — Rtztz = T 4 — ^g-j (3.19) 

Rtxrx R-tyry Rtzrz T ~l~ ^4^2 (3.20) 



Rxyxy Rxzxz Ryzyz ^ ~^ ^4 ' (3-21 ) 



The bulk stress tensor is entirely determined by the contributions to the Riemann 
tensor which do not solve the fluid equations of motion, as it is these that do not solve the 
vacuum Einstein equations. If v is a constant, this yields the boosted black brane, which 
satisfies the vacuum Einstein equations. If v is linear, then again this is a solution of 
the linear order fluid equations as we have checked above, and therefore as we will check 
below v' will not contribute to the gravitational stress tensor at order 0{v). Therefore 
the first nontrivial contributions to the stress tensor arise from the Riemann tensor at 
linear order in v" (v — v' — 0) 

Rtrty = -l(l + -^»i(br)v" (3.22) 

Rxrxy = ^ (2 + 2t>i(6r) + 2&V g'^br) - brX>[(br)) (3.23) 
Rzrzy = \ (2t»i(6r) - br^(br)) (3.24) 
Ryxtx = -4^2(2 + (b 5 r 5 - br)t,[(br)) (3.25) 



i W = jK--"i|- ( 3 - 26 ) 
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As a check on our calculation, we will also calculate the contributions to the various 
tensors at linear order in the nondifferentiated velocity v 

Rtxyx = Rtzyz = £J ^1 ~ V (3.27) 

Rtytr = -( r2 + ^) V ( 3 - 28 ) 

Rtryr = (l ~ V (3.29) 

Rxyxr = Rzyzr = ^ + ^4^2 ^ V (3.30) 

and in v' 

Rtxt y = (b 3 r 3 + &V + (&V - l)(2 fll (6r) + KW))^ (3.31) 

iW = (-2 + &V 4 + (6 5 r 5 + 6r)(2 fll (6r) + brg'^br)))-^ (3.32) 

Rtyxr = -(2 + 6 4 r 4 + (6 5 r 5 + 6r)(2 fll (6r) + brg'^br)))-^ (3.33) 
2t>' 

iW = -(ftV + (6V - l)(2 fll (6r) + 6rfli(6r)))^ (3.35) 

iW = -^ 2 rV(2 'i(^) + &rftf(6r)) . (3.36) 

The Ricci tensor is now easily calculated. Again the order O(v ) terms are those of 
the static black brane solution 

Ru = 4r 2 - ±. (3.37) 

Rt r = R r t — 4 (3.38) 

Rxx = Ryy = Rzz = 4r . (3.39) 

Contributions to the stress tensor will arise from the v" terms, (v — v' — 0) 

Rt y = Ryt = -^(2 + 4(1 + 6 4 r 4 ) Dl (6r) + (6 5 r 5 - br)(t>[(br) + brt>'[{br))) (3.40) 
v" 

R ry = R yr = ^(2 + 4th + br{2br g[(br) - t>[(br) + brX>'[(br)) . (3.41) 
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The terms in the Ricci tensor proportional to v are those of a rigidly boosted black 
brane 

< = ^ *S? = ^ (3-42) 

which provides an exact solution both to the hydrodynamic equations and to Einstein's 
equations with a negative cosmological constant. Again the terms linear in v 1 yield a 
solution to the fluid equations and so Einstein's equations, although only to linear order 
0(v) 

R xy = -(86 3 r 3 01 (6r) + (5&V 4 - l^&r))^- . (3.43) 
Using the large r asymptotic expansions of Ref. [5] 

°. ~ -lab + g|3 + ■ ■ • < 3 - 45 ) 

we find that the asymptotic behaviors of the v" terms in the Ricci tensor are 

13 v" 

~ top? 

^ ~ ■ (3-47) 

The Ricci scalar is 

R = -20 . (3.48) 

There is no contribution at order 0(v) to the Ricci scalar. This is guaranteed for any 
solution of the vacuum Einstein equations with cosmological constant A = —6, and so 
there could not have been any corrections from the v and v' terms. There are no corrections 
from the v" terms at linear order because the corresponding components of the inverse 
metric are themselves of order 0(v), and so the contributions to the Ricci tensor are of 
order 0{y 2 ). 



3.5 Contributions to the Ricci tensor at 0(v ^) and 0(v^) 

Before continuing with the calculation of the bulk stress tensor, we will pause to discuss 
some of the approximations that we have made. We have made two truncations. First, 
we have calculated everything at order 0(v). As we are working in units in which c = 1, 
v is small for nonrelativistic speeds and so this is a valid approximation in a region in 
which the flow is sufficiently slow. 
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A more dangerous truncation is that of higher derivatives of the velocity. The grav- 
ity/hydrodynamics correspondence is a one to one map between gravitational and fluid 
solutions in a derivative expansion. More precisely, the kth order map relates the trun- 
cation of the fluid equations to k derivatives and that of the gravity equations to (k + 1) 
derivatives. The iterative procedure described in Ref. jl] in principle determines this map 
for all k, however in practice this map has only been determined to order k = 2. In other 
words, it provides a metric as a function of v, v' and v ", however a perfect matching with 
Einstein's equations would require also corrections involving the higher derivatives v w 
which are not known. 

General arguments based on dimensional analysis suggest that these corrections be- 
come smaller at higher k. In general one expects that each derivative leads to a contribu- 
tion which is subdominant by a factor of Tl with respect the previous derivative, where 
I is the distance scale of the derivative. Ideally one would like to check this claim for all 
terms with, say, three or four derivatives. However this would require a knowledge of the 
map at orders k = 3 and k = 4. 

The map at order k = 2, which we have used, does produce some terms in the curvature 
which depend on the third and fourth derivatives of v. In this subsection we will verify 
that two of these have the expected convergence scaling, and determine the corresponding 
condition on our fluid flow. In other words, we determine a necessary condition for the 
derivative expansion to apply to our flow. 

The Ricci tensor components R xy and Rt y have corrections from the third and fourth 
derivatives of the velocity respectively 

r{3) = yW(x) (ferp^fer) + Di(6r)) 

^ = _^M (350) 

We want to determine the condition under which R^ is subdominant to K£y . As the 
higher derivatives of v define an interpolating function between two solutions over an 
interval of length d, each derivative is larger than the previous one by about 1/d. In other 
words, d x ~ 1/d. 

To test the sub dominance of , it is sufficient to compare it to the similar term in 

(3) 

R X y , which contains t>i. The ratio of these terms is 

tV rsj ( —L rsj (3.51) 

Rx 3 y rt>( 3 )(x) rd 
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therefore the fourth order term is subdominant if d ^> 1/r in the entire bulk. The bulk 
extends from the horizon at r = 1/6 = nT to the boundary at r = oo. Therefore 
convergence requires 

(3.52) 

This fourth order term is suppressed by irdT with respect to the third order term, in 
line with the above expectations from dimensional analysis. This means that the gravity 
duality procedure is only convergent when d is sufficiently large. Of course, the duality 
never yields a solution of the vacuum Einstein equations, and so one may argue that its 
convergence is immaterial. Nonetheless, it is only well-defined as a series when d satisfies 
f l3~52j) . 



3.6 The static black brane solution 

As a check on our calculation and conventions, we recover that the static (v = 0) black 
brane satisfies the vacuum Einstein equations with cosmological constant A = —6 



R 
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(3.53) 



4 Two Calculations of the Stress Tensor 



In this section we will calculate the bulk stress tensor of the surface layer interpolating 
between the vacuum gravity solutions using two different methods, corresponding to two 
different metrics. First, we will apply the duality map of Ref. [1] to a fluid flow which 
interpolates between the two solutions, the stationary solution on the left and the linear 
velocity solution on the right. In this case, as we have seen, the interpolating region is 
necessarily larger than the inverse temperature. Next, we will directly interpolate between 
the gravitational solutions using the Israel matching conditions |1 1 j . This method requires 
the interpolating region to be very thin, and uses the assumption that in this limit the 
extrinsic curvature remains bounded. 
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4.1 Interpolating between the hydrodynamic flows 



The duality map of Ref. [I] takes a fluid flow and yields a dual metric. This dual metric 
solves the vacuum Einstein equations when the fluid flow satisfies the hydrodynamic 
equations of motion (12.31) . If the flow does not satisfy the equations of motion, the dual 
metric does not satisfy the vacuum Einstein equations. Thus apparently there is no benefit 
in using this map over any other map. However we will use the map, and observe the 
consequences. The resulting dual metric will necessarily solve Einstein's equations with 
some value of the stress tensor 

8-kGnT^ = R^ v - -Rg^ v + Ag^ . (4.1) 

We will determine this value. 

We saw in Eq. (I3.53j) that there is no contribution to the stress tensor at order O(v ). 
We have argued that, at order 0(v), the dominant contributions to the stress tensor are 
proportional to v" . These are easily found from (14. ip to be 

v"{x) (4 (6 4 r 4 - 1) t)i(6r) - br (6 4 r 4 - 1) (t>i(6r) + brt)'{(br)) - 2) 
ty ~ 32nG N b^ ( } 



v"(x) (40i(6r) + br (2brg' l (br) - v[(br) - brv"(br)) + 2 

32-nG N r 2 

There appears to also be a contribution proportional to v' 

v'jx) {br ((6 4 r 4 - 1) g?(6r) + 3br) + (56 4 r 4 - 1) g[{br)) 
xy ~ 16icG N V*r ' ( ' 

At order v' one expects no contributions to the stress tensor, as a solution with a linear 
velocity satisfies the fluid equations at order 0(v). Therefore a nontrivial contribution 
would be in contradiction with the gravity/hydrodynamics correspondence. We will see 
shorty that this contribution is in fact equal to zero. 

The functions t>i(r) and fli(r) are defined as 

2 r .... s r , v-i 



r°° x 3 -1 

0i (r)= / dx——— . (4.6) 



x(x 4 — 1 

Integrating we [5] obtain analytical expressions for t>i(r) and Qi(r 

01 = "I + 5 + 8^ (r4 " :) ('° g £tT* + 2 tan-'(r) - * ) (4.7) 

, I = l(lcg( ( 1 + ^ 1+ ^ -2t m -(r) + ,|. (4. f 
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The derivatives of these expressions are 



0i = o(^— r--^-r + ^-r-r) (4-9) 



r 1 1 

2V^ 2 + 1 r 2 + l r + 1 r 

i, r 2 r 111 

01 = "( r 2 + i)2 + ( r 2 + i) 2 + 2(r2 + 1) + ^ " 2(r + l)2 (4 ' 10) 



and 



t>i = -— (yrr 4 + 2 (r 4 + l) log(r + 1) - 2 (r 4 + l) tan" 1 ) 



-4r 3 + 2r 2 - (r 4 + l) log (r 2 + l) + vr) (4.11) 



t»? = j7 ^t-^ r (3 log (r 2 + 1) + vr(r + 1) (r 2 + l) (r 4 - 3) 

+ fr(r + 1) (r 4 + r 2 -3) -3)r(21og(r + 1) 
- log (r 2 + 1)) - 2(r + 1) (r 2 + l) (r 4 - 3) tan : (r) 

-2 (2r 4 + r 3 + r 2 + r + 3) r 2 - 6 log(r + 1)) . (4.12) 

The explicit formula Eqs. (14.91) and (14.101) for the derivatives of 0i can be combined 
to show that 

r ((r 4 - 1) g'((r) + 3r) + (5r 4 - l) g[(r) = . (4.13) 
This combination is proportional to formula (14. 4 p for T xy , therefore 

T xy = (4.14) 

and there are no contributions proportional to v'. 

Similarly one may evaluate the combination of functions that appears in T ry 

r (2rg[(r) - rX>'[{r) - t>i(r)) + 4t>i(r) + 2 = 0. (4.15) 

This implies that 

T ry = (4.16) 

leaving only T ty , the momentum in the y direction. Thus the bulk stress tensor contains 
no stress, only momentum. 
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We may use the exact expressions for the functions gi and Oi to simplify the only 
nonvanishing component of the stress tensor 

v"(x) 

T ty = - „ l !, 3 . 4.17 

Using the fundamental theorem of calculus, this may be integrated over the interpolating 
region to obtain 

v' 



dx T tv = ~ i« n us 3 ( 4 - 18 ) 

where v' is the derivative of the velocity in the region x > d. In particular, at this leading 
order the integrated stress tensor of the surface layer is independent of the interpolation 
and independent of d. Of course it still depends on the map that we used to generate the 
dual metric. 

Had the v' term been the dominant contribution, the stress tensor would have been 
constant, and so the integral would be have proportional to d. Similarly a v ® term would 
have led to a stress tensor proportional to 1/d, and higher powers of v to other scalings. 
Therefore it is somewhat nontrivial that the leading contribution to the integrated stress 
tensor is in fact d- independent. Clearly this ^-independence is desirable, as d is not a 
physical quantity but merely an artifact of the scheme that we used to regularize the 
divergent second derivative of the fluid velocity. 

The bulk stress tensor does not satisfy the null energy condition. As the only nonva- 
nishing component is T ty , the only nonvanishing product of a null vector w and the stress 
tensor is 

w L Tw = 2w t T ty w y . (4.19) 

As T ty is already of order 0(v), at order 0(v) one need only consider the terms in w of 
order O(v ). That is to say, w only needs to be null with respect to the static black brane 
metric. Consider for example the null vectors w± 



«4 = r, <4 = ±n/l-^j. (4.20) 



The product Q)| is 



v"(x) 

wiTw± = =f — t^t^a/ 1 - T^-r (4.21) 

which is nonzero. However w + and u>_ yield opposite signs, as incidentally do the two 
choices of signs of v. Therefore at least one of these will yield a negative product, and so 
the bulk stress tensor does not satisfy the null energy condition. This may or may not 
mean that no external matter may be consistently added which produces such a surface 
layer. 
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4.2 Israel's matching conditions on the gravity duals 



We will now calculate the bulk stress tensor in a different geometry. Following Ref . [TT] , we 
will consider the vacuum Einstein solution corresponding to a static fluid on the left and 
that corresponding to a linear velocity flow on the right. These solutions will be glued 
together by interpolating continuously between the two metrics over a distance d and 
taking the limit d — > such that the extrinsic curvature remains bounded. In Ref. [TT], 
Israel has shown that the resulting configuration contains two solutions separated by a 
surface layer whose bulk stress tensor is independent of the interpolation used. 

Following Ref. [TT], the first step in the calculation of the stress tensor is the definition 
of the unit normal vector to the hyperplane 

n M = {0,r, 0,0,0}, (4.22) 

which satisfies the normalization condition 

n^ v n v = ^{n x f = 1 . (4.23) 

The surface layer £ extends along all of the directions except for the x direction. A basis 
of tangent vectors to £ is 

ds = e^dx 1 (4.24) 

where 

e (i) = {1,0,0,0,0} (4.25) 

e (tf) = {0,0, 1,0,0} (4.26) 

e w = {0,0, 0,1,0} (4.27) 

e (r) = {0,0, 0,0,1}. (4.28) 

In terms of these tangent vectors the extrinsic curvature may be calculated as 

JCy = e w • V,n = ^ - n m r m>ji = ^ - n m T™ . (4.29) 

On the left, where the fluid is static (v = 0), substituting (I3.16P into ( I4.29P one finds 
no extrinsic curvature 

Ktf = -rl% = (4.30) 
Kj~) = -rT x yr = . (4.31) 
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On the right, where the fluid velocity is linear, the Christoffel symbols of Eq. (13 . 1 T[) yield 
a nontrivial extrinsic curvature K^ + \ 



< ] = = ^3 ( 4 - 32 ) 



Kg> = -rT* yr = V - (1 + b 2 r%(br)) . (4.33) 

The tensor 7^ is defined to be the difference between the extrinsic curvatures on the 
two sides of the surface layer 

7 « = Kip - . (4.34) 
The bulk stress tensor integrated over x is equal to the tensor Sij, defined by 

- torGnSij = i l3 - gi ff™ . (4.35) 

The expression (14.351) for the integrated bulk stress tensor was derived in pj] for a 
4-dimensional space with no cosmological constant. While several factors in the deriva- 
tion change in our current 5-dimensional situation, Eq. (I4.35P remains unchanged. The 
cosmological constant term yields a contribution proportional to the integral of A times 
the metric integrated over the thickness d of the surface layer. As the metric is taken to 
be finite, this term vanishes in the d — > limit. 

The trace of 7 is 0(v 2 ), therefore (I4.35P yields the integrated bulk stress tensor 

S - = -16^( 1 + 6Vfl ''W)' (4 ' 37 > 

These are equal to the integrals over the x directior£]of the stress tensors of Subsec. I4TT1 
at order k = 1, in other words, without the Oi term that entered into the metric (13.21) 
multiplied by v" 



v" 



T »' 1) = -16^( 1 + 6V28 '( fr '')) ' (4 ' 39 > 

The Oi terms arose from the dualization of the interpolating region, which did not satisfy 
the equations of motion. It therefore cannot enter into the Israel calculation, which uses 



1 Note that, following Ref. [TT], the measure of this integral must be that of x rescaled to normal 
coordinates. Therefore the integral contains an additional factor of r = Jlf^ 
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only the solutions of the vacuum Einstein equations. Indeed, the Di terms in (13 .2p are 
singular in the limit d — > as v" diverges as 1/d, and therefore the boundedness of the 
extrinsic curvature assumed in Israel's derivation fails for the metric interpolation (13.21) . 

Like the stress tensor (14.171) calculated by interpolating the hydrodynamic flow, the 
Israel stress tensor does not satisfy the null energy conditions. Again, to linear order in 
v, one may consider vectors which are null with respect to the static black brane metric. 
Therefore, again one may consider the null vectors w± of Eq. (14.201) . As T ty is, at least 
for any finite d, equal to that of Subsec. 14.11 divided by the positive combination br, the 
sign of the inner product (I4.2ip is unchanged. Therefore the null energy condition is also 
violated by this stress tensor. 

The main difference between the two stress tensors is then that T yr does not vanish 
for the Israel tensor. Remembering that in our Gaussian null coordinates the r direction 
is the sum of a spatial and temporal piece, the spatial component implies that there is a 
nonzero stress. More precisely, while both Israel's thin surface layer and the thick fluid 
surface layer have a nonvanishing y momentum, the Israel surface layer also has a flux of 
this y momentum in the radial direction, from the boundary into the horizon of the black 
brane. As the black brane is infinite in the x direction, this is not problematic for the 
time-independence of the solution. 

5 Future directions 

Turbulence often arises as a result of the boundary conditions placed on a fluid. As 
a preliminary step towards an understanding of turbulence in gravity, we have proposed 
two gravitational duals of such boundaries. Both of these duals involve the addition of 
a surface layer of matter, with a certain stress tensor. These proposals are in a sense 
trivial, as the dynamics of the duals is defined not by any known equations of motion, 
but by the duality map itself. It remains to be shown whether such matter can exist. For 
example, even if the equations of motion which it obeys can be found, the existence of a 
UV completion of the matter theory may be fundamentally obstructed as in Ref. [22] . Or 
the failure of the null energy condition may imply that, whatever the ultraviolet theory 
may be, the wall simply disintegrates before it has any significant effect on the fluid. 

Of course, an ultraviolet completion is not necessarily a prerequisite for learning some- 
thing interesting about whatever the gravitational dual to turbulence may be. After all, 
no ultraviolet completion of Einstein gravity is used in this correspondence. The surface 
layer implies the existence of equations of motion which are distinct from the Einstein 
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vacuum equations and perhaps pathological. However the interesting part of the fluid, 
the turbulent part, is not at the wall. For example, if we consider the motion of a fluid in 
a pipe, the flow may be turbulent throughout the interior of the pipe. The ultralocality of 
the duality map implies that, at a distance greater than 1/T from the pipe, the vacuum 
Einstein equations are still satisfied by the gravity dual. Thus in a sense the ultralocal- 
ity decouples the problem of understanding turbulence in gravity from the problem of 
defining a gravity dual of a boundary. 

Besides trying to characterize the gravitational dual of turbulent flow, the other inter- 
esting question is to find the gravitational dual of the conditions under which turbulence 
can occur. In nonrelativistic, incompressible flows, turbulence is expected when the prod- 
uct of a system's characteristic scale L times the characteristic velocity v of a fluid is 
much greater than the kinematic viscosity. In Ref. [8], the authors claim that for the 
conformal fluids dual to AdS black branes, turbulence is expected when LTv ^> 1, where 
T is the temperature of the fluid. The AdS /hydrodynamics correspondence is expected to 
be reliable at scales L such that LT ^> 1. Therefore since v < 1, it appears that whenever 
turbulence is expected, LT > LTv ^> 1 and so the correspondence can be trusted at least 
for quantities that vary over a distance L. (3 + l)-dimensional turbulence is characterized 
by vortices of various sizes from L down to the dissipation scale [10]. Thus the duality 
appears to be reliable at least for the largest vortices in a turbulent flow. The dissipation 
scale is a function of L, T and v, and so in principle one may determine whether or not 
the duality is reliable for vortices all of the way down to this scale and so for the entire 
flow. 

Understanding the gravity duals of turbulent flows, as described above, may yield new 
insights into the dynamics of black branes in AdS space, perhaps revealing a surprising 
difference between branes in AdS4 and AdSs, or indicating that generically they come with 
funnels attached as in Refs. |23J. The main weakness of this program is the dependence on 
asymptotically AdS geometry in the duality map of Ref. [I]. There was no such restriction 
in the original correspondence of Ref. [2] , nor in other identifications of black holes and 
viscous fluids such as the blackfold program of Refs. |23J [25] and the Wilsonian identifi- 
cation of Ref. [2S] ■ An extension of turbulence to asymptotically Minkowski space could 
relate (3+l)-dimensional fluid dynamics to wealth of studies of asymptotically Minkowski 
5d black objects, such as Refs. [27]. More importantly, relaxing the asymptotically AdS 
condition may mean that fluid mechanics, perhaps in only 2+1 dimensions, has something 
to teach us about real world gravity. 
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